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Abstract 

A technique  is  described  for  designing  FIR  band- 
stop  digital  filters  which  effectively  attenuate 
wideband  clutter  in  a digital  Mil  radar  processor. 
The  clutter  return  can  be  caused  by  rain  or  chaff 
and  is  modeled  with  a Gaussian  power-density 
spectrum.  The  bandstop  filter  is  designed  to 
minimize  the  transition  region  for  fixed  filter 
length,  specified  stopband  width,  stopband  atten- 
uation and  passband  ripple. 


This  research  was  supported  by  Army  Research 
Office  Grant  No.  DAAG29-76-G-0212.  . 


Background 

The  primary  Interference  to  a radar  signal  return 
is  called  clutter.  This  clutter  is  typically 
much  stronger  than  the  narrowband  doppler-shif ted 
target  return  and  may  be  caused  by  ground  re- 
flections, rainstorms  or  deliberate  Interference 
such  as  chaff.  The  ground  clutter  is  effectively 
removed  by  a simple  two  or  three  tap  digital 
filter  in  the  moving-target  indicator  (MTI)  radar 
signal  processor.  The  degree  to  which  the  signal/ 
clutter  ratio  is  enhanced  by  the  MTI  filter  is 
called  the  Improvement  factor  (I).  The  standard 
definition  of  I assumes  that  the  signal  gain  is 
computed  by  averaging  over  the  entire  filter 
response  from  dc  to  the  pulse  repetition  frequency 
(1/T)  where  T is  the  pulse  repetition  or  filter 
sampling  Interval.  In  point  of  fact,  the  frequen- 
cy response  of  the  two  or  three  pulse  canceller 
(TPC)  is  quite  dependent  upon  frequency  and  pro- 
vides a true  gain  over  a restricted  portion  of  the 
region  and  severe  loss  over  much  of  the  remaining 
response.  However,  the  average  signal  gain  is 
close  to  0 dB;  hence,  I is  really  a measure  of 
clutter  loss  (C)  which  can  be  more  than  SO  dB  for 
a TPC  and  ground  clutter.  A study  of  the  tradeoff 
between  ground  clutter  rejection  and  usable  pass- 
band  with  more  sophisticated  filters  (N  > 3)  has 
been  reported  by  Houts  & Burlage  [1]. 

Unfortunately,  wideband  clutter  rejection  requires 
a more  sophisticated  filter  oeslgn  because  the 
average  velocity  (v^)  of  the  clutter  spectrum  is 

not  zero  and  has  a standard  deviation  (o^)  which 

is  much  greater  than  the  ground  clutter  value. 
Typically  (2],  for  ground  clutter  < 0.25  m/s 


with  V • 0,  whereas  for  rain  1.0  < a < 3.5  m/s 
c V 

and  chaff  has  1.0  < a <2  m/s  with  v <40  m/s. 

V c 

The  considerable  spread  in  is  due  to  two  facts. 
First,  is  primarily  composed  of  two  independent 
components,  a turbulence  component  and  a 

shear  component  (o  ) which  are  sum-squared  to 
2 ® 

form  ; second,  is  a function  of  range, 
whereas  o^  is  = 1.0  m/s  Independent  of  range. 
Translated  to  C-band  radar  (f^  ■ 5.5  GHz)  this 

implies  the  clutter  spectrum  has  a center 
frequency  (f^,)  between  dc  and  1.5  kHz  and  a 

standard  deviation  (o)  between  37  Hz  and  128  Hz. 

In  the  presentations  which  follow,  the  frequency 
response  will  be  normalized  with  respect  to  the 
sampling  rate  (1/T).  Since  finite-impulse 
response  (FIR)  digital  filters  have  an  amplitude 
response  |H(f)|  which  is  sytmetrlc  with  respect 
to  1/2T,  a normalized  frequency,  F • fT  ranging 
between  0 and  0.5  is  sufficient  to  completely 
describe  the  filter  response. 

FIR  Filter  Design 

The  optimum,  in  a Chebyshev  sense  (minimized 
maximum  error)  , FIR  digital  filter  algorithm  of 
McClellan  et . al.  [3]  serves  as  the  basis  for 
designing  a bandstop  digital  Mil  filter.  Intu- 
itively, the  ideal  bandstop  filter  would  be  de- 
signed by  specifying  the  desired  stopband  center 
frequency  (fg)>  bandwidth  (B^)  and  attenuation 

(A^)  required  to  achieve  the  desired  amount  of 

clutter  rejection  (C) . A practical  filter  must 
also  consider  the  number  of  multipliers  or  taps 
(N) , passband  ripple  (^p)  And  provide  two 

transition  regions  (B^)  between  the  stopband  and 

passband  intervals.  Unfortunately,  the  afore- 
mentioned design  algorithm  hereafter  referred  to 
as  MPAR  utilizes  a relative  error  weighting 
parameter  (W)  in  lieu  of  specifying  values  for  R 

and  A^.  Consequently,  the  design  of  a bandstop 

filter  which  achieves  the  desired  clutter  loss  is 
largely  a matter  of  trial  and  error  as  demon- 
strated in  Table  1, 
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TABLE  1 


MPAR 
(N  - 15,  B^ 

DESIGN 
T = 0. 

1 

RESULTS 

06,  f T - 
s 

0.10) 

RUN 

B T 

W 

R 

A 

t 

P 

s 

1 

0.03 

1 

3.1 

15.0 

2 

0.03 

10 

7.2 

28.1 

3 

0.03 

20 

7.8 

33.5 

4 

0.05 

20 

2.5 

43.0 

5 

0.04 

20 

4.5 

38.0 

6 

0.04 

10 

4.2 

32.5 

The  other  filter  design  parameters  are  assigned 
the  following  values:  N = 23,  B T “ 0.06, 


25  dB  which 
5 dB 


23,  B T 
s 

f^T  =0.1  with  a design  goal  of  C 

Is  to  be  achieved  using  A = 30  dB  and  R 
s P 

centered  at  0 dB  (unit  gain).  The  value  of  R Is 

P 

selected  to  provide  adequate  target  detectability 
throughout  the  passband.  It  Is  apparent  that  ad- 
justments In  Bj.  or  W cause  changes  In  both  R^  and 

A^  and  require  several  trials  before  achieving 

acceptable  performance.  Rablner  [4]  has  reported 
algorithms  for  Iterating  one  particular  MPAR  de- 
sign parameter  until  a desired  error  In  the  pass- 
band  (D  ) or  stopband  Is  achieved.  These 

errors  can  be  used  to  define  A and  R by 

s p ^ 


-20  log  (Dg) 


(1) 


and 


20  log  [(1  + Dp) /(I  - Dp)] 


(2) 


Although  originally  designed  for  lowpass  filters, 
Houts  and  Burlage  [1]  successfully  adapted  the 
technique  for  the  design  of  hlghpass  filters  to 
suppress  ground  clutter  using  B^  as  the  variable. 

The  present  extension  of  this  philosophy  to  band- 
stop  or  bandpass  filters  Is  straightforward  [5]. 
Because  of  the  Inclusion  of  a null  at  dc  or  1/2T 
for  even  values  of  N,  the  application  to  bandstop 
filters  Is  only  practical  for  odd  values  of  N and 
positive  coefficient  symmetry  [6].  Although  In 
the  work  which  follows  N Is  held  constant.  It 
should  be  recognized  that  Increasing  values  of  N 

provide  larger  A or  smaller  values  of  B^  or  R , 
s t p 

all  desirable  features  for  a filter.  The  response 
for  a typical  23-tap  bandstop  filter  designed 
using  the  aforementioned  set  of  design  parameters 
Is  shown  In  Figure  1 along  with  an  Ideal  approxi- 
mation to  this  filter  which  has  uniform  stopband 
attenuation  and  linear  transition  band  response. 


Clutter  Rejection-Idealized  Filter 

A desired  clutter  rejection  can  be  accomplished 
using  some  combination  of  the  parameters  A^ , B^ 

and  B^  assuming  the  stopband  Is  located  at  the 

clutter  center,  i.e. , f • f . What  la  not  clear 
sc 

Is  the  best  choice  for  each  parameter.  The  trade- 
off Is  readily  determined  for  the  Idealized  filter 
shown  in  Figure  1.  The  effect  of  Increasing 


Figure  1.  Frequency  Response  for  FIR  Bandstop 
Filter. 

B^  for  B^  = 1.5o  Is  shown  in  Figure  2 for  various 

choices  of  A . It  follows  that  C = 25  dB  can't 
s 


Figure  2.  Effect  of  Stopband  Attenuation  on 

C-B  Tradeoff. 

8 

be  achieved  for  A^  < 25  dB  and  that  it  requires 

B = 4o  for  A - 30  dE  vs.  3o  for  100  dB.  Of 
s s 

course  a real  filter  design  would  require  a much 
larger  value  of  N for  100  dB  attenuation  than  It 
would  for  30  dB  In  order  to  retain  B^  - 1.5o. 

The  effect  of  changing  B^  is  shown  In  Figure  3 

for  A ~ 30  dB.  Again,  It  follows  that  C • 25  dB 
s 

can  be  achieved  either  with  B^  - 2(j,  • 3.1o  or 
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Figure  3.  Effect  of  Transition  Bandwidth  on  C-B 
Tradeoff.  ® 

9 

3^  = lo,  B = 4.5a,  hcwever,  the  total  unusable 
t s 

bandwidth  (B)  defined  by 

B = 2[B  + 2B^]  (3) 

s t 

Is  14.2a  In  the  first  case  and  only  13.0a  In  the 
second  case.  It  Is  also  evident  that  the  wider 
transition  widths  (small  N)  play  a significant 
role  In  clutter  loss  when  the  stopband  Is  narrow 
(B^  < 4a)  and  have  a negligible  effect  when 

B > 6a.  The  effect  on  clutter  loss  (C)  of  off- 
s 

setting  the  clutter  from  the  center  of  the  stop- 

band  (f  ^ f ) Is  shown  In  Figure  4.  The  nonnal- 
c s 

Ized  frequency  offset  Is  defined  as 

Af  »|f  - f|/f  X 100%.  (4) 

c ' s c'  s 

A ten  percent  frequency  offset  reduces  C by  3 dB 
while  a 20  percent  offset  costs  8 dB.  Since  It 
would  be  Impractical  to  continuously  match  the 
filter  to  the  changing  average  clutter  velocity 

(v  ) It  follows  that  the  filter  must  be  designed 
c 

with  more  clutter  rejection  than  needed  under 
Ideal  conditions  to  protect  against  a change  In 
v^,  e.g. , a 10  dB  margin  permits  a 25  percent 

change  In  v^.  Due  to  the  symmetry  of  the  filter 

and  clutter  spectrum  the  direction  of  velocity 
change  Is  Immaterial. 

Clutter  Rejection  - FIR  Filter 

A series  of  15  and  23-tap  FIR  bandstop  filters 
' were  constructed  using  the  transition  bandwidth 
(B^)  as  a variable  parameter  adjusted  to  meet  the 

specifications  A^  • 30  dB  and  R^  ■ 5 dB.  A 15-tap 


Figure  4.  Effect  of  Frequency  Offset  on  Clutter 
Loss . 

filter  was  designed  with  f^T  ■ 0.10  and  clutter 

rejection  data  Is  presented  In  Table  2 assuming 
Gaussian  clutter  with  aT  • 0.01  or  0.02.  The 
Improvement  Is  relatively  constant  for  small 
B^T  (<0.04)  which  can  be  explained  by  noting  that 

Bj.  Is  decreasing  with  Increasing  B^  in  such  a 

manner  Chat  the  unusable  bandwidth  (B) , defined 
by  Equation  (3),  remains  roughly  constant. 
Furthermore,  little  additional  clutter  rejection 
Is  possible  once  B^  > 4a. 

TABLE  2 

15-TAP  FILTER  RESULTS 


(A  - 30  dB» 
s 

R » 5 dB,  f T 
P s 

- 0.10) 

B T 

9 

B^T 

C(oT  - 0.01) 

C(oT  - 0.02) 

0.00 

0.052 

22.3 

13.6 

0.01 

0.048 

23.3 

13.8 

0.02 

0.043 

25.6 

14.4 

0.04 

0.068 

27.1 

21.5 

Data  for  clutter  rejection  of  a 23-tap  filter  is 
presented  In  Table  3 for  f T • 0.15  and 

aT  • 0.02  or  0.03.  Again  C Is  less  for  the 
larger  value  of  o.  What  might  appear  surprising 
at  first  Is  to  note  that  with  aT  • 0.02,  C is 
larger  for  a 15-tap  filter  than  for  a comparable 
23-tap  design.  The  explanation  Is  that  B^  is 

roughly  50%  larger  for  the  15-tap  design,  hence 
a larger  portion  of  the  clutter  spectrum 
experiences  some  attenuation.  Naturally  this  Is 
achieved  at  Che  expense  of  additional  unusable 
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b andwld  th . 


TABLE  3 

23-TAP  FILTER  RESULTS 


(A^  - 30 

dB,  R - 5 dB, 
P 

f T - 0.15) 
s 

B T 
s 

B T 
t 

C(aT  » 0.02) 

C(aT  - 0.03) 

0.00 

0.030 

6.9 

3.7 

0.02 

0.029 

10.3 

5.8 

0.04 

0.040 

18.5 

11.1 

0.06 

0.036 

22.4 

13.2 

0.08 

0.031 

26.4 

15.0 

The  tradeoff  between  B and  C la  demonstrated  In 
Figure  5 which  clearly  shows  that  a particular 
level  of  clutter  rejection  requires  larger  B for 
broader  clutter  spectra.  The  23-tap  data  adheres 
closely  to  the  theoretical  results  obtained  from 

an  Ideal  filter  h = 1.5o  and  variable  B . 

t s 

This  Is  somewhat  expected  due  to  the  similarity  of 
the  frequency  responses  shewn  in  Figure  1.  It  is 
evident  that  B^  > 1.5a  for  the  15-tap  data  pre- 
sented In  Table  2.  Closer  adherence  between  the 
Ideal  and  actual  filter  clutter  rejection  could  be 
obtained  by  varying  N until  • 1.5o,  a result 

which  was  closely  approximated  by  much  of  the  23- 
tap  data  shown  In  Table  3.  Essentially  Identical 


c 


Figure  5.  Clutter-Unusable  Bandwidth  Dependence  on 
OT. 

B^  results  were  obtained  for  a 23-tap  filter  de- 
signed with  f^T  • 0.1.  The  effect  of  offsetting 
the  clutter  center  frequency  (f^)  from  the  center 
of  the  stopband  was  studied  for  two  different 


23-tap  filters  (B^  = 1.5o  and  B^  = 4o).  For  each 

study  the  clutter-loss  data  were  equivalent  to 
the  theoretical  results  shown  in  Figure  4. 

Conclusions 

It  is  possible  to  achieve  reasonable  wideband 
clutter  rejection  with  an  FIR  bandstop  digital 
filter  which  Is  optimal  In  the  Chebyshev  sense. 

An  Idealized  filter  can  be  used  to  examine  the 
clutter-unusable  bandwidth  tradeoff  and  select 
reasonable  design  values  for  A^,  B^  and  B^.  If 

desired,  the  value  of  B^  can  be  adjusted  by 

changing  the  number  of  taps.  The  amount  of  clut- 
ter rejection  using  23-tap  filters  compares 
favorably  with  results  for  the  idealized  filter 
using  similar  stopband  characteristics.  Increased 
clutter  rejection  can  be  achieved  by  increasing 
the  stopband  width  or  stopband  attenuation  thus 
requiring  larger  N to  maintain  constant  B. 
Additional  Investigations  are  required  to  select 
the  center  frequencies  and  numbers  of  bandstop 
filters  required  to  effectively  cover  the  wide 
range  of  possible  clutter  velocities.  Although 
the  filter  design  algorithm  is  not  tied  to  the 
Gaussian  power-density  assumption,  the  idealized 
filter  clutter-bandwidth  tradeoff  curves  from 
which  the  filter  design  parameters  are  selected  ' 
are  dependent  upon  this  assumption  and  would  re- 
quire recomputation  for  another  clutter  spectrum. 

) 
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